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Abstract 

£NJ ' We consider corrections to nucleon correlation functions arising from times that are far from the 

asymptotic limit. For such times, the single nucleon state is contaminated by the pion-nucleon 
and pion-delta continuum. We use heavy baryon chiral perturbation theory to derive the spectral 
representation of the nucleon two-point function. Finite time corrections to the axial current 
correlation function are also derived. Pion-nucleon excited state contributions drive the axial 
correlator upward, while contributions from the interference of pion-delta and pion-nucleon states 
drive the axial correlator downward. Our results can be compared qualitatively to optimized 
nucleon correlators calculated in lattice QCD, because the chiral corrections characterize only low- 
energy excitations above the ground state. We show that improved nucleon operators can lead to 
D ' an underestimation of the nucleon axial charge. 
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I. INTRODUCTION 



Lattice gauge theory simulations continue to make impressive progress towards addressing 
quantitatively the non-perturabative regime of QCD [l|. Two decades ago, it was thought 
impossible that lattice QCD simulations would confront experimental data without orders 
of magnitude of increased computing power, and considerable algorithmic advances. Today, 
however, state-of-the-art simulations have pushed forward on all fronts by: including fully 
dynamical quarks, shrinking the lattice spacing, enlarging the lattice volume, and decreasing 
the size of the light quark masses (including a sprint to the physical point [2]). The current 
status of these simulations is reviewed in |3j. Lattice QCD is approaching a point of not 
only complementing existing experimental programs, but potentially guiding future ones. A 
thorough overview of hadron structure from recent and forthcoming lattice calculations is 
presented in 0]. 

Performing lattice QCD calculations with light pions is an exciting recent development, 
but one that is accompanied by new problems. At fixed lattice sizes, for example, finite 
size effects from virtual pion fluctuations will become increasingly important. Another 
notorious issue in dealing with hadrons on the lattice is the signal-to-noise problem. In 
the chiral regime, this problem will become rather acute. Consider the nucleon two-point 
correlation function, G(t). Over long times, r, the two-point function has an exponential 
falloff governed by the nucleon mass, M^. Over such long times, however, the statistical 
noise in the correlator, S(r), is dominated by its coupling to three-pions [5|. This leads to 
the signal-to-noise ratio having the behavior 



with iV as the number of independent measurements on gauge configurations. As pion masses 
enter the chiral regime, the long-time behavior of the correlation function is dominated by 
statistical noise. Generally this limits lattice QCD measurements of hadronic correlation 
functions to times that are not ideally long. 1 

In this work, we consider corrections to hadronic correlation functions arising from the 
pion continuum. The Lehmann-Symanzik-Zimmerman (LSZ) reduction formula provides the 
field theoretic recipe for producing hadronic states in the limit of asymptotically long times. 
Over shorter times, however, there are transient fluctuations to multiparticle states within 
a hadron. As the signal-to-noise problem restricts lattice QCD simulations to times far 
from the asymptotic regime where LSZ applies, we are motivated to investigate corrections 
to the reduction formula. To focus our discussion, we consider the chiral dynamics of 
the nucleon. We derive corrections to nucleon two- and three-point functions using chiral 
perturbation theory. In heavy baryon chiral perturbation theory, one expands about small 
energies above the nucleon mass. The nucleon operator should thus be thought of as some 
baryon interpolating field that has been highly optimized through the removal of high- 
lying modes, i.e. modes k that are k > A x above the nucleon mass, where A x ~ 1 GeV 
is the chiral symmetry breaking scale. In the language of lattice QCD, the renormalized 
nucleon operator in chiral perturabtion theory behaves similarly to a quark-smeared baryon 



Excluding pion zero modes by using parity-orbifold boundary conditions has been suggested as a means 
to overcome the signal-to- noise problem [y, |7| . 
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interpolating field. It is obviously difficult to compare our results quantitatively with those 
from similarly optimized nucleon interpolating fields currently used in lattice QCD studies. 
Our results, however, provide a clean field-theoretic example of the difficulties encountered 
from excited state contamination. Additionally qualitative information, such as the sign of 
corrections from excited states, may be garnered from our computations. For three-point 
functions, the sign is non-trivial due to the lack of a spectral representation. We find that 
contributions from pion-nucleon states are driving the axial correlator up, while those from 
the interference of pion-delta and pion-nucleon states drive the correlator down. 

The organization of our presentation is as follows. First in Section [Til we derive the 
spectral representation of the nucleon propagator in heavy baryon chiral perturbation theory. 
Here we investigate the effects from pion fluctuations on the time dependence of nucleon 
two-point functions. Next in Section IHfl we extend our consideration to the LSZ reduction 
of the nucleon three-point function involving the axial- vector current. Sizable corrections to 
the axial correlator are found, and can lead to overestimation of the nucleon axial charge, Ga- 
On the other hand, we describe a scenario in which the axial charge can be underestimated 
by improving the nucleon overlap in two-point functions. Concluding remarks are given in 
Section M\ 



II. SPECTRAL REPRESENTATION 

Consider a free fermion field, N(x), of mass M. In the limit of large mass, the velocity 
becomes a good quantum number in coordinate space: and, to leading order in 1/M, we 
can write down a heavy fermion effective Lagrangian B B E3, C = NW 4 N, where N is a 
two-component Pauli spinor. The coordinate space propagator, D(x,0), for the free heavy 
fermion has a simple form, D(x, 0) = 6(x4)5(x)e~ Mxi . When interaction terms are included 
in the Lagrangian, the full two-point function, G(x,0), has the general form 



G{x,0) 
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6{x A )5{x)e- Mxi / dEp{E)e~ Ex \ (1) 
Jo 



which follows from taking the large mass limit of the Kallen-Lehmann spectral representa- 
tion. The spectral function, p(E), is a positive function in the distribution sense. Assuming 
there is an isolated single-particle state corresponding to the N, we write p(E) = 5(E)+p(E). 
The residual spectral function, ~p(E), is then assumed to vanish below the energy E t ^, with 
-E t h > 0. With the single-particle state isolated, the spectral representation has the form 



G{x,0) = 6{x 4 )8{x)e 



-Mxa 
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/ dEp{E)e~ Exi 

JE th 



(2) 



In the absence of bound states, E th is the threshold energy to create a multiparticle state. 

Suppose there is an excited state, N*, of mass M* contributing to the spectral function. 
The spectral weight around the value E = M* — M then contains a delta-function. For 
asymptotically large time separation between source and sink, x± 1, this excited-state 
contribution is suppressed relative to the ground state by an exponentially small factor: 
exp[— (M* — M)x 4 ]. The contributions from multiparticle states too are suppressed in the 
limit of asymptotic time separation. If we assume the spectral weight near threshold is of 
the form ~p{E) oc (E — i^)" -1 , then the contribution from the multiparticle branch-cut is 
suppressed by a factor of (x 4 )~ n exp(— E th Xi), for large time separation. The power n is 
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FIG. 1: One-loop diagrams contributing to the nucleon two-point function. Single (double) lines 
denote nucleons (deltas), while the dashed lines denote pions. The filled circles denote axial 
couplings from the interaction Lagrangian, Eq. ([3|). 

determined by the available phase space at threshold, and the required angular momentum 
for multiparticle production. If correlation functions are not deduced at large enough time 
separations, the light pions will likely present difficulty for bound states in QCD, because 



A. Chiral Computation 

Treated as heavy fermions, free nucleons can be described the heavy fermion Lagrangian, 
with N now upgraded to a two-component isospinor, iV = (p, n) T . Interactions between 
pions and nucleons; as well as pions, nucleons and deltas can be included in a way consistent 



with chiral symmetry [11|, [12j. These interactions, moreover, can be systematically organized 
in terms of a small expansion parameter e, where the quantities: k/M with k as residual 
baryon momentum, p/A x with p as pion momentum and A x = Airf as the chiral symmetry 
breaking scale, m n /A x , and A/A x with A as the nucleon-delta mass splitting, are all treated 
to be of the size e. The leading interactions among the baryons and pions are at 0(e), and 
contained in the Lagrangian 

C = 2g A N^S ■ AN + g^ N [T f • A N + N ] A ■ T] . (3) 

Here S = cr/2 is the spin operator, and the quartet of deltas are packaged in a Rarita- 
Schwinger field T. The axial-vector field of pions is at leading order A = j + . . ., with 
= T a ir a . 

Using Eq. ([31), we can determine the interacting nucleon two-point function at one- loop 
order in the chiral expansion and thereby determine the spectral function. The contributing 
diagrams are shown in Figure [TJ Their computation is complicated by divergences which we 
regulate using dimensional regularization. To renormalize the two-point function, we choose 
to work in coordinate space. 2 This is particularly advantageous for our consideration of the 



2 It is instructive to carry out the computation additionally in momentum space. After the mass and 
wavefunction renormalization have been taken into account, the renormalized self-energy W en (E) has the 
form 

£ ren (£)=£(£)-E(0)-i?^£(0), 

but still contains divergences. These divergences only contribute to the extreme short-time behavior of 
the two-point function, i.e. contributions to G(r) in Eq. (j4|) of the form S'(t), and 8"(t). The spectral 
function, however, can be deduced from the regulated self-energy from analyticity. The relation 

P(E) = --Qtm [G(E)}, 
produces the same spectral function as deduced in coordinate space. 
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LSZ reduction for three-point functions which we undertake in Section IHII 
We first define the zero-momentum projected two-point function 



GO 



dx G(x, r; 0, 0), 



(4) 



where we have used spacetime translational invariance to locate the nucleon source at the 
origin of our coordinate system. To renormalize the zero-momentum two-point function, we 
require that there be a single nucleon contribution at asymptotically large r. This consists 
of two parts, the mass renormalization 
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(5) 



that fixes the nucleon mass to its physical value (which thus absorbs its pion mass depen- 
dence), and the wavefunction renormalization 



lim [ 



(6) 



that fixes the single nucleon probability to unity. Carrying out the one-loop chiral compu- 
tation of the nucleon two-point function, we find 3 



G(r) 
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where the pion-nucleon and pion-delta fluctuations are described the by spectral weights 



m 



213/2 



(4tt/) 2 E 2 
l6gl N [(E - A) 2 - m 2 ] 3 / 2 



E 2 



(8) 
(9) 



For asymptotically large times, the decaying exponential from the single nucleon state is 
the dominant contribution to the two-point function. Modification from pion interactions can 
be deduced in this limit by considering the spectral weights near threshold. The pion-nucleon 
and pion-delta spectral functions both vanish as the 3/2-power of the energy available at 
threshold, cx (E - E th f' 2 . This power-law behavior is due to the two-body phase space, cx 



3 It is straightforward to derive analogous results in a finite spatial volume. With periodic boundary 
conditions, the lattice momentum modes are quantized in the form k — 2iTn/L, where n £ Z 3 , and L is 
size of the lattice, which is assumed to be the same in each of the three spatial directions. Integrals over 
the energy are replaced by sums over the momentum modes 

dEEf(E) - £ ^/(v/fcHmJ), 

n 

where n = \fr\? . We are assuming a large enough volume, m^L 3> 1, so that the sum can be replaced by 
the integral for which the spectrum is continuous. Finite volume corrections to this approximation can 
be calculated as the difference between the sum and integral. 
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FIG. 2: Interacting nucleon effective mass plots. The curve denoted by Full corresponds to the 
effective mass derived from the fully interacting two-point function in Eq. (|7|), while that denoted by 
No A is the same two-point function evaluated without delta resonance contributions, i.e. gatv = 0. 
The vertical bar symbolizes the region where we imagine noise might dominate the correlation 
function. 



(E — EtYi) 1 ^ 2 1 and the requirement that the pion and baryon be in a relative p-wave. Because 
the pion-delta threshold is larger than the pion- nucleon threshold by the mass splitting A, 
the pion-delta contributions are suppressed relative to pion-nucleon contributions by the 
exponential factor exp(— Ar). Hence we shall neglect the delta contribution in asymptopia. 

Expanding the pion-nucleon spectral weight about threshold, we arrive at the asymptotic 
expansion of the nucleon two-point function 



G(r) 



phys 7 



1 + V27T 



4tt/ J (m^r) 5 /2 



1 
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8 m^T 128 (m n Ty 



+ 



(10) 



The . . . denotes power-law suppressed terms proportional to l/m w r. The expansion for large 
times is only an asymptotic one, which is evidenced by the particularly large numerical 
coefficients of higher-order terms. The times over which the expansion in l/m n T results 
in a controlled approximation to Eq. ([7]) are quite large, r > 50 a, for a = 0.125 fm. For 
smaller times, the first term in the expansion yields the best agreement with Eq. (|7j), as is 
characteristic of asymptotic expansions used outside their range of validity [13]. To analyze 
nucleon two-point functions, however, we shall not use the approximation in Eq. ffTUl) . but 
return to the full form in Eq. (]7j). 



B. Mass Extraction 

Now we consider the effect of the pion continuum on the nucleon two-point function. A 
typical quantity used for guiding the eye in fits of lattice two-point functions is the effective 
mass, and we shall employ this tool to investigate the interacting nucleon two-point function 
derived from chiral perturbation theory. The effective mass function is defined as 

M eff (r/a) = - log G{T n + a \ (11) 
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We employ lattice units using a temporal lattice spacing of a = 0.125 fm, which is the spacing 
in some current-day lattices. For large r/a, the effective mass should become flat with a 
value corresponding to the nucleon mass. The onset of a plateau in the effective mass gives 
one an indication of when the ground state dominates the two-point correlation function. 
Of course, in actual lattice calculations the signal-to-noise problem for the nucleon limits 
one to times that are not ideally long. 

In Figure [2J we make two effective mass plots using the interacting nucleon two-point 
function in Eq. (j7]). The low-energy constants are fixed to their known values, specifically we 
use g A = 1.27, g A N = 1.5, A = 290 MeV, and / = 130 MeV. We plot the ratio of the effective 
mass to the nucleon mass for two values of the pion mass, m n = 140, and 300 MeV. The 
onset of a plateau in the effective mass hence corresponds the ratio M^/Mn approaching 
unity which is exhibited in both panels of the figure. We see, moreover, that the pion- 
delta continuum gives negligible contributions due to the exponential suppression factor, 
exp(— At). In a recent lattice study [14], the noise begins to dominate the signal around 
r/a = 12 at m n ~ 300 MeV. We show a vertical line at this time to denote an imagined 
noise barrier for our chiral computation. Fitting the nucleon two-point function with a 
single exponential form at = 300 MeV for 6 < r/a < 12 will overestimate the nucleon 
mass by ~ 3%. For the plot at the physical pion mass, we assume that the imagined noise 
barrier remains at r/a = 12. Fitting the two-point function on a smaller time interval, say 
8 < r/a < 12 results in the same quantitative overestimation of the nucleon mass. We can 
see, moreover, that improving the signal-to-noise at light pion masses will expose effective 
masses which continue to fall with time. These qualitative features can be deduced directly 
from the positivity of the spectral representation. 

One final observation about the nucleon two-point function is as follows. Improving the 
overlap with the nucleon in two-point functions largely removes contamination from the 
pion-nucleon continuum, as evidenced by Figure [2j Due to the exponential suppression, 
pion-delta contributions are largely unaffected by the improvement. Such contributions 
are not necessarily suppressed, however, in other nucleon correlation functions. The axial 
correlator provides such an example. 



III. AXIAL THREE-POINT FUNCTION 

We now treat the nucleon axial-current correlation function at finite times using the 
coordinate-space approach. Formally there are no spectral representations for three-point 
functions. Consequently contamination from excited states can drive the correlator up or 
down depending on the underlying dynamics. 



A. Chiral Computation 



Nucleon three-point functions are constructed using the LSZ reduction formula to iso- 
late single-nucleon contributions from the external legs. We are interested in a three-point 
function formed from inserting the axial current J^(y,i) between two nucleon states. Fol- 
lowing the typical lattice procedure (for recent lattice calculations of the axial charge, 



sec 



15l . |16|, [l7j, [18|, [19J), and projecting both source and sink onto vanishing three mo- 
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FIG. 3: One-loop diagrams contributing to the axial-vector current matrix element of the nucleon. 
The open squares denote axial-vector current insertions, while the remaining diagram elements are 
as in Figure [TJ 



mentum, we form the ratio 



Jdxjdy (0\N(x, r)J+(y, *)JV+(0, 0)|0) 



U ^ T ^ fdx(0\N(x,T)M(O,0)\0) ' (12) 

Using the axial current at tree-level in the effective theory, namely = 2g A N^S^t+N, we 
arrive at 

K 5li (r, t) = 6{t - t)6(t) 2g A u^S.u, (13) 

where u is a Pauli spinor. Working beyond tree level, there are divergences. After regular- 
ization, we renormalize to the physical axial charge, G A , namely 

n 5tx (r,t) 2G A u^S,u, (14) 

which additionally absorbs the pion-mass dependence of the axial charge into the physical 
coupling G A . 

The leading one-loop diagrams contributing to the three-point function of the axial cur- 
rent are shown in Figure [3J The required terms of the baryon axial current appear as 

J+ = 2g A iV^r+iV + g AN (T^r+iV + iVV+Tj - 2g AA T^S fl r + • T, (15) 

while the pion-nucleon and pion-nucleon-delta interaction terms are contained in the La- 
grangian Eq. fl3]). Additionally to calculate the ratio TZ 5 ^(r,t), we must divide by the 
two-point function calculated to one-loop order. The result of the chiral computation can 
be cast in the form 

TZ 5 ,(r,t) = 2u^S,u [G A + F A (T,tj\e(r-t)9(t), (16) 

where the function F a (t, t) encodes the deviation from the limit of infinite time separation. 
The form of this function is 

F A (r,t) = ^ f°°dEp Nn (E) [e- E ^) + e -« - e -^] 

-^f £dE vffiU^[e- £M (1 + e-*<-«>) 

+ e~ Et (1 + e- At ) - e~ Er (e' At + e~ A ^) ' 

+ (g A + f dE-p^E) [e~ E ^ + e- E * - e^] . (17) 

The nucleon-pion, p N7T (E), and delta-pion, p A7T (E), spectral weights have been given above 
in Eqs. (JSJ) and (jUJ). In the limit {r, t} — > oo with r > t, accordingly we have F A (r,t) — > 0. 
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B. Axial Charge Extraction 



To investigate the effect of finite times on the extraction of the nucleon axial charge, we 
plot the ratio of axial three-point to two-point functions, R 5fl (r,t), given simply by 

R 5fl (r,t) = G A + F A (r,t). (18) 

Scaling R 5 ^(r,t) by the axial charge, Ga, we thus expect a plateau at unity for asymptot- 
ically large time separations, r ^> t ^> 1. We fix the sink time at rja = 12, which is a 
typical source-sink separation in lattice QCD computations of three-point functions. With 
the sink time r fixed, the current insertion time dependence of is shown in Figure H] 
We compare the behavior with and without delta resonances. Including the resonances re- 
quires knowledge of the delta axial charge, #aa- This low-energy constant is poorly known, 
see (2Q| | . and we assume the value #aa = —2.25. Rather fortuitously the time-dependence 
of the axial correlator is largely insensitive to the value of the delta axial charge. To a 
good approximation, we can drop all terms in Eq. (ITTj) that arise from diagrams with only 
intermediate state deltas, that is 

F A (r,t) « ^ rdEp m {E) [e- E ^+e- Et -e^] 

-^f [ dE vasmj^H (1 + e -A M) 

+ e~ Et (1 + e~ At ) - e~ Er (e" Ai + e~ A ^) ] , (19) 

which is independent of #aa- The figure confirms that Eq. ([19]) well approximates the time- 
dependence of the axial correlation function. Excluding the delta completely (by setting 
9 an — 0), however, has a dramatic effect on the time-dependence of the axial correlation 
function. The figure shows that the contributions from pion-nucleon excited states drive the 
axial correlator upward. Those from pion-delta states are negligible, while contributions from 
the interference of pion-delta and pion-nucleon states drive the axial correlator downward. 

We see the effect of insufficient time can lead to spurious plateaus for three-point func- 
tions. At both pion masses depicted, the derived correlation function flattens out, but 
above the asymptotic value of unity. Our expectation from chiral perturbation theory is 
that the nucleon axial charge would be overestimated, with the overestimation worsening as 
the pion mass is lowered. Rather large couplings of the pion to the nucleon, and nucleon- 
delta transition complicate extraction of quantities from lattice three-point functions. As 
such systematic errors, moreover, depend on the pion mass, the resulting chiral behavior of 
Ga, for example, will be specious. Lattice QCD calculations often extract the ratio of the 
axial-to- vector charges, Ga/Gv- The vector charge, Gy, is determined from the three-point 
function of the vector current which is also susceptible to an effect from insufficient time. 
The one-loop calculation of the vector three-point function in heavy baryon chiral perturba- 
tion theory, however, does not produce any finite time corrections due to exact cancellation 
borne in by the vector Ward identity. This cancellation is only guaranteed at zero momen- 
tum transfer. Thus away from this limit, the vector three-point function too is subject to 
finite-time corrections. 

As a final comment, there has been debate about underestimation of the axial charge in 



lattice computations [24], |25| . The most recent study of the axial charge in fully dynamical 
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FIG. 4: Time-dependence of the axial current correlation function. For sink time r/a = 12, we 
plot the ratio of the axial three-point to two-point correlators as a function of the current insertion 
time t/a. This ratio we scale by the axial charge G a- The curve denoted by Full corresponds to 
the axial correlation function determined from Eq. (|17|) . while No A corresponds to the correlator 
without delta-resonance contributions, i.e. gAN = 0. The curve No ExpA excludes from the axial 
correlator contributions that are suppressed by an exponential factor involving the mass splitting 
A, see Eq. (fT9l) . 



three-flavor simulations [18( finds Ga falling as the pion mass is decreased. This behavior 



was attributed to a finite volume effect, but not that described by effective field theory [21|, 



22, 23 



4 We undertook the present study with the hope that the falling value of Ga could 
be attributed to excited state contamination: i.e. although the axial charge is calculated 
in 
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on two volumes, these volumes have differing temporal extents, and hence differing 
effects from insufficient source-sink separation. Given our results, there is a plausible way 
in which such underestimation is possible. Underestimation of the axial correlation function 
can arise from improvement of the nucleon overlap in two-point functions. Designing baryon 
operators with the best ground-state overlap in a two-point function does not necessarily 
optimally improve three-point functions. In a two-point function, improvement of the the 
nucleon operator decreases contamination from the pion-nucleon continuum, but not from 
the the pion-delta continuum, because the latter makes only negligible contributions. In the 
axial three-point function, however, the dominate excited state contributions arise from both 
the square of the overlap with pion-nucleon states, and the interference term between pion- 
nucleon and pion-delta states. Improving the nucleon operator diminishes the former faster 
than the latter. As the latter is negative, the axial correlator can thus be underestimated 
for highly optimized nucleon interpolating fields. 



IV. SUMMARY 

In an interacting quantum field theory, single particle states emerge only in the limit of 
long times. The existence of particle interactions necessarily implies that at shorter times a 



4 



While the volume dependence employed in [18| is inspired by the p-regime of chiral perturbation theory. 



the size of the corrections is not consistent with expectations 21 1 . 
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single particle state can fluctuate into multiparticle configurations. Such configurations can 
affect the extraction of hadronic observables. In particular, the existence of light pions in 
the spectrum of QCD complicates the extraction of hadronic properties as they introduce a 
characteristically long time scale, r ~ l/m w , over which pion branch-cuts can be important. 
Using the nucleon in chiral perturbation theory as an example, we investigated multiparticle 
contributions to two- and three-point functions. Specifically we found the nucleon mass is 
subject to overestimation, which is directly related to the positivity of the spectral function. 
For the nucleon axial charge, Ga, we found that three-point functions calculated away from 
asymptotic time separations are affected by the competition between pion-nucleon and pion- 
delta contributions. Pion-nulceon excited states drive the axial correlation function upward, 
while the interference between pion-delta and pion nucleon excited states drives the axial 
correlation function downward. The sign of these corrections depends on the underlying 
chiral dynamics, and cannot be deduced from general principles alone. 

The formulae derived in this work apply to the renormalized nucleon operator in chiral 
perturbation theory. Because high- lying excitations are absent, this operator behaves like a 
quark-smeared lattice QCD interpolating field. While not quantitative, the analytic insight 
offered by our computation, however, is useful for comparing qualitatively with lattice QCD 
computations, nearly all of which use optimized baryon operators with smeared quark fields. 
Using a basis of interpolating operators, a variational method can be employed to isolate 



ground states from contamination by excited states [26|, |27|. Work along these lines has 
been pursued, for example, the baryon spectrum has been investigated in jH, H^, [3(J 31 . 
Removing pion-nucleon contamination from the nucleon two-point function should improve 
the mass determination, although the study in 30j seems to overestimate the nucleon mass. 
Despite eliminating a large part of the excited-state contamination, the nucleon effective 
mass continues to drop (their Figure 6 has the same qualitative behavior as our Figure [2]), 
signaling the need for a two-exponential fit. Only a single exponential fit, however, was 
carried out. 

Our calculation of the nucleon axial current makes plausible that the use of improved 
nucleon operators can lead to an underestimation of the axial charge. The extraction of 
nucleon observables from three-point functions is potentially fallible due to the pion contin- 
uum, which is important owing to the lightness of the pion, as well as the comparatively 
large axial coupling, g&. The same observation can be made for three-point functions in- 
volving deltas due to the size of the axial transition coupling, g^N- As pion masses approach 
the physical point, the effect of insufficient time in lattice QCD calculations of nucleon and 
delta properties will become pronounced. Optimizing the overlap with the ground state in 
two-point functions will not necessarily lead to better determination of observables from 
three-point functions because of interference terms. We hope our work can qualitatively 
guide the extraction of nucleon and delta observables from lattice QCD. 
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